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We define quantum cellular automata as infinite quantum lattice systems with discrete time dy- 
namics, such that the time step commutes with lattice translations and has strictly finite propagation 
speed. In contrast to earlier definitions this allows us to give an explicit characterization of all local 
rules generating such automata. The same local rules also generate the global time step for au- 
tomata with periodic boundary conditions. Our main structure theorem asserts that any quantum 
cellular automaton is structurally reversible, i.e., that it can be obtained by applying two blockwise 
unitary operations in a generalized Margolus partitioning scheme. This implies that, in contrast to 
the classical case, the inverse of a nearest neighbor quantum cellular automaton is again a nearest 
neighbor automaton. 

We present several construction methods for quantum cellular automata, based on unitaries com- 
muting with their translates, on the quantization of (arbitrary) reversible classical cellular automata, 
on quantum circuits, and on Clifford transformations with respect to a description of the single cells 
by finite Weyl systems. Moreover, we indicate how quantum random walks can be considered as 
special cases of cellular automata, namely by restricting a quantum lattice gas automaton with local 
particle number conservation to the single particle sector. 
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I. INTRODUCTION 

The idea of generalizing the classical notion of cellular 
automata to the quantum regime is certainly not new. 
Indeed, it is already present in Feynman's famous paper 
from 1982, in which he argues that quantum com- 
putation might more powerful than classical. However, 
although there have been several formal definitions of 
quantum cellular automata over the years, the theory is 
not in good shape at the moment, and a systematic ex- 
ploration of the general properties of such systems on the 
one hand and of the potential for computational applica- 
tions on the other has hardly begun. We believe that this 
is partly due to deficiencies of the existing approaches, 
and therefore propose a new one, which is very natural, 
and only requires a few basic assumptions: a discrete cell 
structure with a finite quantum system for every cell and 
translation symmetry, a discrete time step for the global 
system, reversibility, and finite propagation speed. 

Quantum cellular automata ("QCAs") are of inter- 
est to several fields. There are obvious connections to 
the statistical mechanics of lattice systems, and poten- 
tial applications to ultraviolet regularization of quantum 
field theories. In Quantum Computer Science they ap- 
pear as one natural model of computation extending the 
well-developed theory of classical cellular automata into 
the quantum domain. But also the experimental side 
is rapidly developing: quantum computing in optical 
lattices j2[ and arrays of microtraps [3( are among the 
most promising candidates for the first quantum com- 
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puter that does useful computations £|. It is typical for 
such systems that the addressing of individual cells is 
much harder than a change of external parameters af- 
fecting all cells in the same way But this is just the 
theoretical description of a cellular automaton. Possi- 
ble tasks, which have the same built-in translation in- 
variance are simulations of solid state models. Typically 
classical simulations run into problems already for mod- 
erate systems sizes, precisely because of the dimension 
and complexity explosion which Feynman noted, and for 
which he proposed quantum computation as a cure. The 
theory presented in this paper can be considered as pro- 
viding the first elements of an assembly language for such 
simulations. 

The problems which have plagued previous attempts 
to define QCAs begin with the definition of a system of 
infinitely many cells. Consider the simplest operation 
such an automaton should be able to perform: applying 
the same unitary transformation separately to each cell. 
This would involve multiplying infinitely many phases, so 
there is really no well-defined unitary operator describing 
the global state change. Therefore, the quantization ap- 
proach "just make the transition function unitary" does 
not work very well. Similarly, the notion of state vectors 
as amplitude assignments to uncountably many classical 
configurations is ill-defined. But this would be a can- 
didate for the "configurations" of a QCA, which causes 
problems for a definition of QCAs in terms of configu- 
rations and their transformations. Various approaches 
HI EJ IE 13 will be described and commented in Sec- 
tion [V] A constructive method to obtain QCAs, which 
is common to most of these approaches (including ours) 
is partitioning the system into blocks of cells, applying 
blockwise unitary transformations, and possibly iterating 
such operations. Model studies based on such construe- 
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tions (e.g., 01) therefore produce results independently 
of the definition problems. 

In order to arrive at a satisfactory notion of QCAs, it is 
helpful to draw on ideas from a discipline, which has been 
dealing with infinite arrays of simple quantum systems 
for a long time, i.e., the statistical mechanics of quantum 
spin systems. Infinite systems have been considered par- 
ticularly in the algebraic approach to such systems [ll| . 
The basic idea is to focus on the observables rather than 
the states, i.e., to work in the Heisenberg picture rather 
than the Schrodinger picture [12J. The main advantage 
is that in contrast to localized states, it does make sense 
speak of local observables 0], i.e., observables requiring 
a measurement only of a finite collection of cells. The 
global transition rule of a cellular automaton is then a 
transformation T on the observable algebra of the infi- 
nite system. As always in the Heisenberg picture, the 
interpretation of such a transformation is that 'prepar- 
ing a state, running the automaton for one step, and 
then measuring the observable A' gives exactly the same 
expectations as preparing the same state and measuring 
T(A). As always, T must be completely positive and sat- 
isfy T(l) = 1. But more importantly we can state the 
crucial localization property of QCAs: When A is local- 
ized on a region A of the lattice, then T(A) should be 
localized in A + N = {x + n\ x e A, n 6 N}, where M is 
the neighborhood scheme of the QCA. 

To our knowledge, this view of QCAs was first used in 
[l4| . where the approach to equilibrium in a QCA with 
irreversible local rules (based on a partitioning scheme) 
was investigated. The general picture of QCAs remained 
unsatisfactory, however, because the partitioning scheme 
seemed a rather special way of constructing a QCA. For a 
satisfactory theory of QCAs we demand that there should 
be a direct connection between the global transition rule 
T and the local transition rule: If we know the global 
transition rule, we should be able to extract immedi- 
ately the local rule in a unique way, and conversely, from 
the local rule we should be able to synthesize the global 
rule. The class of global rules should have an axiomatic 
specification, the most important of which would be the 
existence of a finite neighborhood scheme in the above 
sense. On the other hand, for the local transition rules 
we would prefer a constructive characterization. That is, 
there should be a procedure for obtaining all local rules 
leading to global rules with the specified properties, in 
which all choices are clearly parameterized. 

The partitioning QCAs of [14| failed to meet these re- 
quirements, because they provided a construction, but no 
axiomatic characterization of the global rules obtained 
in this way. In particular, it remained unclear whether 
two steps of such an automaton could be considered as a 
single step of an automaton with enlarged neighborhood 
scheme. The idea enabling the present paper was that 
all these difficulties vanish if we restrict to the class of 
reversible QCAs. The axiomatic characterization is ex- 
tremely simple: In addition to the above locality condi- 
tion we assume that the global rule must have an inverse, 



which is again an admissible quantum channel. This is 
equivalent to saying that T must be an automorphism of 
the observable algebra. Then the local rule is simply the 
restriction of this automorphism to the algebra of a single 
cell. Conversely, since every observable can be obtained 
as a linear combination of products of single-cell observ- 
ables, the local rule determines the global automorphism. 
This allows us to subsume all the known constructions of 
QCAs, but also to prove a general structure theorem: 
every reversible QCA is structurally reversible, i.e., we 
can write the local rule in a partitioning scheme involv- 
ing two unitary matrices, which makes it apparent how 
the global rule can be unitarily implemented on arbitrar- 
ily large regions, and how to obtain the local rule of the 
inverse. 

A further bonus from our proof of the structure the- 
orem is that it does not actually require the global rule 
to be an automorphism: it works under the prima facie 
much weaker assumption that the global rule is a homo- 
morphism (and not necessarily onto). Then invertibility 
follows (see Corollary below). Therefore, invertibility 
was not included in Definition 1, which makes it much 
easier to verify whether a proposed rule is indeed a QCA. 

The structural invertibility was an open problem in the 
theory of classical reversible cellular automata in higher 
dimensional lattices until recently |15| . Hence, since our 
proof of the quantum result is rather simple, it appears 
that some proofs in the classical domain can be simplified 
by going quantum. 

Our paper is organized as follows: In Section [H] we 
begin with the axiomatic definition of QCAs in the sense 
described above. Its counterpart, the constructive de- 
scription is given in the form of a collection of basic con- 
structions and examples in Section IlIII It turns out that 
one of these constructions, based on partitioning, is al- 
ready sufficient to obtain all QCAs in the sense of our def- 
inition. This rather surprising result is stated and proved 
in Section [Tvl The ideas of the proof also allows us to 
give an explicit parameterization of the simplest class of 
QCAs: nearest neighbor automata in one dimension with 
one qubit per cell. As mentioned in the introduction, the 
current literature on the subject is mostly based on a def- 
inition we do not find satisfactory. We discuss these, and 
some further related definitions, in more detail in Sec- 
tion Finally, in an appendix we provide some math- 
ematical background on finite dimensional C*-algebras, 
which play a key role in the proof of Theorem |SJ 



II. DEFINITION OF QCAS 

We consider an infinite cubic array of cells, labelled by 
integer vectors x € Z s , where s > 1 is the spatial di- 
mension of the lattice [0| . Each cell contains a d-level 
quantum system with the same finite d > 2. That is to 
say, with each cell i e Z s , we associate the observable 
algebra A x of the cell, and each of these algebras is an 
isomorphic copy of the algebra of complex d x d-matrices. 
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When A C If is a finite subset, we denote by A(A) the 
algebra of observables belonging to all cells in A, i.e., the 
tensor product ® xeA A x . By tensoring with unit op- 
erators on A 2 \ Ai we consider .4(Ai) as a subalgebra 
of ^4(A2), whenever Ai C Aa< In this way the product 
A1A2 of Ai G A(Ai) becomes a well-defined element of 
A(Ai U A2). Moreover, tensoring with the identity does 
not change the norm, so we get a normed algebra of lo- 
cal observables, whose completion is called the quasi-local 
alaebra\n\. and will be denoted by A(Z S ). Similarly, for 
other infinite subsets A C Z s we define -4(A) as the clo- 
sure of the union of all A(A') with A' C A finite. 

When x G If is a lattice translation, we denote by 
t x the isomorphism from each A y to A x + y , and its ex- 
tensions from .4(A) — » _4(A + x), by shifting every site. 
Here we have used the notation A + .t = {y + x\y G A} for 
shifted lattice subsets, which we also extend to Ai + A 2 = 
{x x + x 2 \xi e AJ. 

A state oj of the spin system is a linear functional on 
A(Z S ), which is positive in the sense that u>(X*X) > 
and normalized as w(l) = 1. Equivalently, a state is 
given by a family ui\ of density operators on (C d )® A (for 
each finite A), such that u(X) = tr(w A A) for X G .4(A). 
The local density matrices have to satisfy the consistency 
condition that, for Ai C A2, WAi is obtained from lo\ 2 
by tracing out all tensor factors in A2 \ Ai. Note that a 
state does not correspond to a configuration of a classical 
automaton, but rather to a probability distribution over 
global configurations. 

Definition 1 A Quantum Cellular Automaton with 
neighborhood scheme Af C Z s is an homomorphism 
T : A{1f) — > A(Z S ) of the quasi-local algebra, which 
commutes with lattice translations, and satisfies the lo- 
cality condition T(A(A)) C A(A + J\f) for every finite set 
A C Z s . The local transition rule of a cellular automa- 
ton is the homomorphism T : Aq — ► A{Af) . 

Note that a unitary operator for the time evolution 
is not necessary in this formulation. Instead we have 
replaced it by its action on observables. Of course, the 
time step has to be read in the Heisenberg picture. That 
is, measuring some local observable A G .4(A) at time 
t + 1 is equivalent to measuring the observable T{A) at 
time t. The transition rule thus describes this backwards 
calculation for a single cell. To(.4o) is some isomorphic 
copy of the one-cell algebra embedded in a possibly quite 
complicated way into the algebra of neighboring cells. 
The relationship between the global evolution and the 
one-site transition rule is as simple as it should be: 

Lemma 2 

(1 ) The global homomorphism T is uniquely determined 
by the local transition rule Tq. 

(2) A homomorphism Tq : Aq — > A(Af) is the transition 
rule of a cellular automaton if and only if for all x G 
If such that J\f n (N + x) ^ the algebras To(.4o) and 
t x (To(Ao)) commute elementwise. 



Proof. By translation invariance the action of T x : 
A x — > A{M + x) is determined as T X (A X ) = t x Tqt_ x (A x ) 
on all other cells. Moreover, because T is a homomor- 
phism, the extension to any local algebra is also fixed. 
Explicitly, consider a product ® xeA A x = JlxeA °f 
one-site operators. This equation just expresses our iden- 
tification of the one site algebras A x with subalgebras of 
^4(A) by tensoring with unit operators. Then the homo- 
morphism property of the global evolution requires that 




Note that the product on the right hand side cannot be 
replaced by a tensor product, because the factors have 
overlapping localization regions x + TV. Moreover, the 
argument of T is a product of commuting factors, hence 
so is the right hand side. Hence the the commutativity 
condition (2) is necessary. 

Conversely, if the factors T X {A X ) commute, their prod- 
uct is unambiguously defined. Since every local observ- 
able can be expressed as a linear combination of tensor 
products, Eq. defines a homomorphism on the local 
algebra, as required. This shows the converse of (2), and 
since we have given an explicit formula of T in terms of 
To it also shows (1). □ 

It is clear that the commutation condition of the 
Lemma can be expressed as a finite set of equations, and 
can therefore be verified effectively. Since only a small 
portion of the lattice is needed in this verification, the 
same steps are needed to check local transition rules fore 
QCAs on graphs which locally look like If . Such graphs 
can be seen as integer lattices with periodic boundary con- 
ditions. The only condition we will need to impose is that 
the periods for the boundary condition are not too small 
compared with the size of neighborhood scheme N (a 
condition called "regularity" below) . No algebraic condi- 
tions on the homomorphism To are needed. Therefore we 
can turn the construction around, and immediately get 
a QCA on the infinite lattice from a QCA with finitely 
many cells. Since for finitely many cells a homomorphism 
is always just implemented by a unitary matrix, this al- 
lows us to define QCAs even for the infinite system, just 
by specifying unitary matrices with suitable properties. 

Let us describe the periodic boundary QCAs more pre- 
cisely and see what conditions are needed for the neigh- 
borhoods. The cells of a system with periodic boundary 
conditions arise from the cells in If by identifying cer- 
tain cells, namely all those differing by a vector 7 in some 
subgroup r C If . The set of cells is thus identified with 
the quotient C = If /T. For each point in x G C, i.e., 
each equivalence class x — xq + T of identified cells, only 
one observable algebra is given. The sites C can be rep- 
resented in a so called fundamental domain of T like the 
parallelogram in Fig. ^ But there seems to be no way to 
draw this nicely as a set of square cells. 

For points x = (xq + T) G C, i.e., we can define the 
translation x + n = xq + n + T . Consider now a neigh- 
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borhood scheme TV C If . The neighborhood of x G C is 
then x + Af C C. Note that as far as the model with pe- 
riodic boundary conditions is concerned we could change 
each n € Af by a lattice vector in T without changing the 
neighborhood of any point. But since we are interested 
in the connection with the infinite model, we do not ad- 
mit this ambiguity. The neighborhood scheme is called 
regular for the given periodic structure given by £, if the 
equations we have to check for the commutation rule of 
a cellular automaton are the same in both cases. Since 
these equations depend only on the intersections between 
translates of TV we only need to make sure that the geom- 
etry of intersections is the same. So suppose that neigh- 
borhoods on C intersect, say (x + Af) n (y + Af) / 0. 
This means that there is a translation m£Z s such that 
x + m — y and also TV n (to + TV) ^ 0. Clearly, the first 
condition determines to up to a lattice translation T, and 
the second is an intersection condition on the infinite lat- 
tice. Obviously, x + (Af n (TV + to)) c (x + TV) n (y + Af) . 
But the inclusion could be strict if there is more than 
one to with the required properties. This is precisely 
the case regularity must exclude. To summarize, we call 
a neighborhood scheme TV C If regular for a subgroup 
r C If , if TV n (to + TV) ^ and n 6 T, n =f 0, imply 
that TV fl (to + n + Af) — 0. A more compact equivalent 
form is (TV + TV - TV - TV) n T = {0}. An example of a 
regular neighborhood is given in Figure^ That the same 
neighborhood becomes non-regular for a smaller lattice 
is shown in Fig. [21 




FIG. 1: Periodic boundary conditions. The parallelogram rep- 
resents the given finite lattice C, where cells cut by opposite 
boundaries must be suitably joined. The marks on the corners 
join up to a full circle. Two translates of the same neigh- 
bourhood are shown. It is regular, because the geometry of 
intersections is the same as on the infinite plane. 

Then checking condition (2) of the Lemma for the QCA 
on C and for the QCA on If are exactly equivalent. We 
call this useful principle the Wrapping Lemma: 

Lemma 3 The QCA transition rules on a finite lattice 
C with respect to a regular neighborhood scheme Af are 
in one-to-one correspondence with the transition rules for 
QCAs on If with the same neighborhood scheme. 
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FIG. 2: The same neighborhood scheme is not regular on a 
5x4 torus: Again we have two translates, but their intersec- 
tion (cross-hatched) cannot be realized on the infinite square 
lattice as intersection of two such neighborhoods. 

III. BASIC CONSTRUCTION METHODS 

A. Commuting Unitaries and Phases 

Consider a unitary operator Uq in some local al- 
gebra A(Af), which commutes with all its translates 
U x = t x (Uq), up to a phase. That is, we require that 
there are complex numbers Cx with = 1 such that 
U T X (U Q ) = ( x t x (Uo)Uo for (Af + x) H TV + 0. This is 
equivalent to 

U X Uy — C,y — X Uy TJ X ( 2 ) 

We can then formally define a unitary operator 

-U= l[U x ". (3) 

Here the the scare quotes indicate that there is no way 
this infinite product can be made sense of. However, we 
can define instead the action of this "operator" on local 
observables. To this end, note that the actions A i— > 
U*AU X commute for different x. Moreover, if x+Af does 
not intersect the localization region of A, this action is 
the identity. Therefore, in the infinite product of these 
operations only a finite set is not the identity, and their 
product defines an automorphism T. More formally, we 
have 

T(A) = Urn U* A AU A , (4) 

A/"Z 

where A £ A(Z S ), and U A = H xeA U x . The limit is 
over any sequence of finite sets, eventually absorbing all 
lattice points, and if A is localized in a finite region, the 
limit is actually constant for sufficiently large A. 

The local transition rule is found by applying T to a 
single cell. This gives the neighborhood scheme 

Af = \J{ti + x\OEAf + x}=Af-Af. (5) 

X 

We mention three special cases of this construction: 

• When TV = TV = {0}, we have to choose a unitary 
operator U a G Aq acting on a single cell. Thus 
the cellular automaton acts by applying the same 
unitary rotation separately to every cell. 
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• Fix a basis in C and consider any unitary Uq 
which is diagonal in the corresponding product ba- 
sis. Clearly, this guarantees that Uq commutes with 
its translates. That is, we can generate a QCA by 
an arbitrary choice of local phases in some compu- 
tational basis. A prominent example is an Ising 
interaction 

H = (1 - 0-3) ® (1 - £T 8 ) (6) 

turned on for a suitable finite time (e.g., t = 7r/4), 
which generates the initial entangled state of a one- 
way quantum computer |l7j together with a cell- 
wise Hadamard rotation. 

• Choose at every site a family of unitaries V p , p — 
1, . . . , d commuting up to a phase. Examples are 
the Pauli matrices for d = 2, or a discrete Weyl 
system. Then Uq can be any finite product of op- 
erators from this family localized on neighboring 
cells. The automata generated in this way form a 
group. 

For all QCAs constructed from commuting unitaries we 
have the strange property that they show no propaga- 
tion, in the sense that the localization region does not 
increase when we iterate the automaton. The reason is 
that n steps are implemented by a product of commut- 
ing unitaries, each of which appears n times. Hence it is 
exactly equivalent to take instead a single step with the 
basic unitary operator Uq replaced by Uq. 

Of course, most QCAs do have propagation. For ex- 
ample, we could take a single step constructed from com- 
muting unitaries, followed by a site-wise rotation along 
a skew axis. Another rich class of propagating QCAs is 
given in the next section. 



B. Quantization of Classical Reversible CAs 

A typical feature of the local phase automata is that 
they leave invariant the algebra of operators diagonal in 
the chosen basis. This algebra V is the quasi-local alge- 
bra of a classical CA embedded into the quantum system, 
which has d states per cell, when A x is the algebra of 
d x d-matrices. It is therefore natural to look for cellular 
automata, which leave this classical subalgebra T> invari- 
ant as a set, but not elementwise. Clearly, such QCAs 
induce a classical CA on the classical subsystem. In fact, 
every reversible classical CA can be obtained in this way. 

Intuitively, this is seen as follows: the classical CA can 
be run with periodic boundary conditions, for simplicity, 
so we have only a finite set C of cells. It then defines a 
permutation of the classical configurations, which we 
can interpret immediately as a unitary permutation oper- 
ator U. This unitary operator is now used to implement 
the local evolution of the QCA. All we need to verify is 
that for A G Aq the evolution U*AU is indeed contained 



in some local algebra A(Af) for some regular neighbor- 
hood J\f. Then the wrapping Lemma asserts that the 
QCA is also well-defined on the infinite lattice. 

However, this argument must be more subtle than it 
looks: it is well known, that the injectivity of a classical 
CA on the infinite lattice implies the existence of an in- 
verse CA but the inverse is hard to compute, because 
there is no a priori upper bound on its neighborhood size. 
Superficially, the inverse neighborhoods do not seem to 
enter the above argument. However, the argument for 
the locality U*AU E A(Af) requires more than the local- 
ity of the classical rule and the unitarity of U. Consider, 
for example, the rule 

c a; +1 = c x + l + c x + c x-l l (7) 

where c* E {0, 1}, and addition is mod2. With peri- 
odic boundary conditions of length L this is an invertible 
transformation unless L is divisible by 3. This proviso is 
not of the form "for sufficiently large L..." , which means 
that the classical automaton does not allow a local inver- 
sion, i.e., there is no inverse cellular automaton. By the 
wrapping Lemma, it is clear that for this rule we can- 
not find a quantum version either [Tflj. This shows that 
the local invertibility of the classical CA must enter the 
argument. We therefore assume now that the classical 
CA is locally invertible, and the lattice C is chosen suffi- 
ciently large, so that the neighborhood schemes for both 
the classical automaton and its inverse are regular for C. 

The classical configurations are functions a : £ — > A, 
where A — {1, . . . , d} denotes the set of classical states for 
each cell, and at the same time labels the computational 
basis of C d . We will write a E A c , and denote by a x = 
a(x) the value of the cell x in configuration a. 

Extending this to a product basis, each configuration 
a E A c determines a basis vector \a). Then the global 
unitary transition operator is defined by U\a) = \F(a)), 
where F denotes the global classical transition function. 
The transition rule of the QCA is determined by com- 
puting all matrix elements of the operator To(|co)(eo|), 
i.e., 

<o|3b(|co>(eo|)|6) - (a\U*(\c )(e \®t c ^)U\b) 

= {F(a)\(\co)(e \®t C \ {a} )\F(b)) 

This expression is = 1, if 

if F(a) = c , F(b) = e , 
and F(a) x = F{b) x for x ^ (8) 

and = otherwise. We have to show that this is of the 
form X <g) l®AA/\ 

Lemma 4 Let F be a classical cellular automaton with 
neighborhood scheme Mc, which has an inverse automa- 
ton with neighborhood scheme A/j. Then T as defined 
above is a QCA with neighborhood scheme J\f — Mc — 
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Of course, this Lemma only gives an upper bound on 
the size of the neighborhood scheme. Depending on the 
particular automaton, A/" may be much smaller. 

Proof: We have to show that for all eo, Co, the operator 
r (|co){e |) is of the form X ® w i tn X E A(M). 

This is equivalent to two conditions: on the one hand the 
matrix elements (a|7o(|co)(eo|)|fr) must vanish, whenever 
a x ^ b x for some point x ^ M, and, moreover, the value 
of the matrix elements must be independent of a x for 
such x. 

Suppose that the matrix element is non-zero, i.e., con- 
dition (JSJ) holds. Then for all y such that y +Mi ^ the 
computation of the the values of a y and b y from F(a) 
and F(b) will give the same values of a y and b y . In other 
words, the diagonality condition holds for y (—A//). 

The dependence of the matrix element on a x is also 
governed by condition JSJl : since the matrix elements can 
only be or 1, we have to show that the validity of the 
condition does not depend on a x for x ^ J\f, given the 
range of equality of a's and b's established in the previous 
paragraph. Indeed, once that diagonality is established, 
one can see that most of the conditions become re- 
dundant. If x is such that [x + Mc) H (—A/}) = 0, then 
F(a) x and F(b) x are computed by the local rule of F 
from identical data, so they must be equal. It there- 
fore suffices to consider the condition for those finitely 
many a; for which a dependence remains possible, i.e., 
x E (— Mc —Mj). But then only those a y enter, which 
contribute to F(a) x via the local rule (and similarly for 
b). This restricts y to (Mc — Mc — Mi), and this is what 
the Lemma claims as the localization region. □ 



C. Partitioning 

The easiest way to build a cellular automaton with 
readily verified locality properties is based on the ccllwisc 
unitary rotations, with the modification of changing the 
partitioning of the system into cells [2(|. The typical 
construction would thus be: 

1. possibly divide the given cells into suitable subcells, 
by writing the one-cell Hilbert space C d as a tensor 
product of other spaces. 

2. partition the set of cells of the previous step into 
blocks in some periodic way: every cell now be- 
longs to exactly one block, and any two blocks are 
connected by a lattice translation. 

3. apply the same unitary operator to each block al- 
gebra. 

4. possibly repeat this procedure with different block 
partitions 

5. possibly split and regroup once again to come back 
to the original pattern of d-dimensional cells. 



Every step in this construction is well-defined for the 
global system for the same reason that cell- wise rotations 
are valid QCA operations. Moreover, the unitary oper- 
ators doing each of the steps are essentially arbitrary, 
and by just inverting the steps we can immediately con- 
struct the inverse QCA. This is why partitioned cellu- 
lar automata are sometimes called structurally reversible. 
Moreover, the partitioning idea allows one to construct 
QCAs from irreversible local rules just as easily as re- 
versible ones 0, 0| . 

Our main theorem (Theorem below) will tell us that 
every QCA can be written in partitioned form. For near- 
est neighbor automata it is sufficient to take two steps 
in which cells are grouped in cubes of side 2, possibly 
with a choice of unequal cell sizes in the intermediate 
step (see Section TlV A() . This is known as the Margolus 
partitioning scheme (see Figure |3J). 



FIG. 3: The Margolus partitioning scheme in s = 2 dimen- 
sions. Operations are alternatingly applied to the solid and to 
the dashed partitioning into 2x2 squares. The square shape 
of the cells is irrelevant, as they only serve to label localized 
quantum systems. 



D. Circuits 

Of course, it is natural to think of a cellular automaton 
as a physical device, which just happens to be infinitely 
extended (or periodically closed). This suggests build- 
ing QCAs from some basic supply of circuit elements, 
whose properties will then ensure that the overall op- 
eration makes sense. The mathematical details of the 
description must then work out automatically, because 
"Hardware cannot lie" . Consider the following example 



FIG. 4: A proposed QCA circuit. 
Here the flow of information is from top to bottom. 
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The symbol stands for a CNOT gate, by which the bit on 
the line with circle and cross (the "target bit" ) is flipped 
if and only if the value of the bit on the line with the fat 
dot (the "control bit") is "1". This is a standard gate 
also for quantum computation. We can readily compute 
the action of this device on classical information: each 
output bit depends only on the input on the same line 
and the line one step to the right. Presumably, this would 
also be the description of the action on computational 
basis states in the quantum case. So can we not just 
build this device, and construct its proper mathematical 
description along with the hardware? 

The problem with this automaton becomes apparent 
already when we try to compute its classical inverse: this 
requires at each switch to know the value of a control bit, 
which is not yet determined. In fact, the inverse does not 
exist, because the initial states "all 1" and "all 0" are 
both mapped to "all 0" . So this device is not reversible, 
even classically. (Incidentally, this also holds for periodic 
boundary conditions, so it is not a problem of infinite 
size). What went wrong? The problem is the timing of 
the gates. In fact, in the usual gate model of quantum 
computation it is assumed that each gate is executed at a 
certain time. Here the times overlap, and if we insist on 
the gates being executed, say from the left to the right, 
we either need infinitely many operation times per step, 
or we run into time ordering problems at the boundary 
condition. 

One way to avoid this problem is to insist on some fi- 
nite number of clock cycles per QCA step, that each gate 
is executed in one of these cycles, and that the gates run- 
ning in the same cycle do not access the same registers. 
An "unscrambled version" of the above impossible QCA 
is drawn in Fig. [S] 



o- 



FIG. 5: An operational QCA circuit taking two clock cycles 

Clearly, this is a partitioning QCA and, in fact, the de- 
scription we just gave of a circuit with timing constraints 
is nothing but the definition of a partitioning QCA. 



E. Clifford automata 

In the implementations of quantum computation there 
is often a separation between "easily implemented" oper- 
ations and others, which may be more costly. For exam- 
ple, "linear" transformations on the quantum light field 
can be performed with mirrors, beam splitters and phase 
plates, whereas squeezing or photon number counting are 
more costly. A similar choice of a subgroup of "easy" op- 
erations for qubit quantum computation is the group of 
transformations, which take tensor products of Pauli ma- 



trices into tensor products of Pauli matrices, the so called 
Clifford arouvffi^ . Indeed in some implementations these 
play a special role, and can be executed in parallel [l7j . 
Clearly, it is important to understand this subgroup com- 
pletely, although it is also clear that such transformations 
alone will not allow quantum computational speedup. 

The natural mathematical setting for the investigation 
of Clifford QCAs are discrete Weyl systems, with a finite 
Weyl system acting at each site, and the tensor products 
of one-site Weyl operators generating a Weyl system with 
infinitely many degrees of freedom. If the local Weyl sys- 
tem has prime dimension d, one can give a very explicit 
description of the group of Clifford QCAs, which will be 
presented elsewhere ..22]. For illustration let us just take 
qubit automata d = 2 in one dimension. 

What is needed to define such an automaton? Since 
the Pauli matrices o~ x and o~ z generate the one site al- 
gebra, we only need to specify the two operators Tq{u x ) 
and To((7y). The Clifford property means that these two 
operators must be tensor products of Pauli matrices, so 
we have 



Tq((t x ) 
To(a x ) 



a(0 (9) 



where each can take the values 0, x,y, z, with 

(Jo = 1- Now the two strings £ = (f-jv, • • • , £iv) and 
rj = (tj_ ft, . . . ,T]ft) completely characterize the automa- 
ton. But which strings are allowed? From the general 
theory we immediately get the necessary and sufficient 
conditions: cr(£) must commute with all its translates, 
the same holds for (7(77). Moreover, <r(£) commutes with 
Ti(a(rf)) for i ^ and anti-commutes for i = 0. Since 
tensor products of Pauli matrices always either commute 
or anti-commute, one can run a computer search for all 
examples with low neighborhood size N. This turns up 
the surprising result that (up to a common translation) 
the strings £ and 77 must be palindromes, i.e., = 

The general theory |22| confirms this, and moreover 
etablishes an isomorphism of the group of Clifford QCAs 
with the group of 2 x 2-matrices, 



(10) 



whose entries are polynomials over the two-element field 
F2 = {0, 1} in one indeterminate z, such that the deter- 
minate is the constant polynomial 1. Here the coefficients 
of £± are bit strings which together determine the string 
£o,£i, used in Eq. ©. 

One can also find a simple set of generators: Apart 
from one-site transformations and the shift only one QCA 
is needed: 



T (a x ) = 1 <g> a z <g> 1 
T (a x ) = a z ®a x ®a z , 



(11) 



possibly, however, acting not between neighboring sites 
as written here, but between sites at a fixed distance L, 
so that the chain breaks up into L non- interacting chains. 
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The prototype has a number of interesting prop- 
erties. For example, the iterates T'(cr(C)) of an initial 
Pauli product <r(£) have a specific form: they consist of a 
Pauli product moving to the left, and another one moving 
to the right, each at maximal speed, and the expanding 
space between these patterns is filled by one of four possi- 
bilities: all 1, all Oy or an alternating patterns oi o~ x ®o~ v , 
or the same shifted by one cell. 

IV. STRUCTURE 

In this section we will employ the commutativity prop- 
erty of transition rules to get some information about the 
structure of possible rules, aiming at the proof that all 
QCA transition rules can be understood in a partitioning 
scheme. 

Without loss of generality, we consider only nearest 
neighbor rules on a cubic lattice. For a non-cubic lat- 
tice we can choose a family of basic cells (a "fundamen- 
tal domain" ) such that all cells are generated from these 
basic ones by translation symmetries. By considering 
the fundamental domain and its translates as new cells, 
we effectively get a family of lattice cells labelled by Z s . 
If the neighborhood scheme involves more than nearest 
neighbors, we can again enlarge cells. Of course, these 
operations partly destroy the underlying lattice symme- 
try, so that operations on regrouped cells may fail to have 
the translation (or other) symmetry of the original lat- 
tice. However, for the proof of structural invertibility a 
regrouped lattice of "supercells" works just as well. 

We begin by describing the geometry of the general- 
ized Margolus partitioning scheme. In the following sub- 
section we state the main theorem: this scheme indeed 
suffices for all QCAs. The proof relies on the concept of 
"support algebras", and is described in Subsection II V CI 
When the support algebras are abelian, one can charac- 
terize the possible QCAs at the single cell level, without 
partitioning (see Subsection IIV E|l . This allows us to de- 
termine explicitly all nearest neighbor qubit automata in 
one dimension fSection flV Gfl . 

A. Generalized Margolus Partitioning 

Consider a cellular automaton with one-site algebra 
A® = Md, lattice Z s , and nearest neighborhood scheme: 

M = {x e Z s | Vil^i | < 1} . (12) 

We will use a cell grouping introduced by Margolus • 
In this scheme one "supercell" is the unit cube 

□ = {x e Z s \ Vi^ = 0,1} , (13) 

which consists of 2 s cells. The even translates D+2x with 
x G Z s cover the whole lattice. We denote by Q the set 
of 2 s quadrant vectors q G Z s for which each component 
is Qi G {— 1, +1} (see Fig[SJ). In particular the vector into 



the positive quadrant (with all Qi = +1) will be denoted 
by 1. Note that the sum of two quadrant vectors is an 
even lattice translation in 2Z S . Moreover, the 2 s cubes 
□ + q for q G Q are disjoint and their union contains all 
neighborhoods of cells in □. Just like the even translates 
of □, the cubes □ + 1 + 2x with x G Z s form a partition 
of the lattice (compare Fig. [2J • 
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FIG. 6: The Margolus scheme in s = 2 dimensions. The 
basic cube "□" is at the center, with the origin marked "0". 
The shaded squares are reached from the origin by quadrant 
vectors. Dashed outlines mark the copies □ + q shifted by 
quadrant vectors. 

A partitioned automaton based on the Margolus 
scheme would alternatingly apply blockwise unitary 
transformations to the cubes in the two partitions. But 
not every QCA can be written in this way. A good coun- 
terexample is the shift in the quadrant direction 1. Here 
the entire quantum information in the Z s cells □ will 
have to be moved to □ + 1 although these blocks have 
only a single cell as overlap. It turns out, however, that 
a slight generalization of the Margolus scheme suffices to 
represent every QCA: all we have to do is to allow dif- 
ferent cell sizes in the intermediate step. For the rest of 
this subsection we will explain the resulting scheme. 

With each quadrant vector q G Q we associate an ob- 
servable algebra B q c A(D + q), which is isomorphic to 
the algebra of n(q) x n(g)-matrices for some integer n(q). 
Since B q is contained in a local algebra, it makes sense 
to consider its (even) translates T x (B q ) C A{D + q + x) 
with x G 2Z S . In particular, A(D + 1) contains all the 
algebras Ti- q (B q ). A crucial assumption of our construc- 
tion is that these subalgebras of A(0 + 1) commute, and 
together span A(D + q). This is possible if and only if 
the equation 

\{n{q) = d^ (14) 

<?6Q 

holds for the matrix dimensions. Note that each cell 
A(D + x) has this dimension, whether or not x is even or 
not. 

The local rule of an automaton now defines (and is 
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defined by) a homomorphism 



(15) 



Since the dimensions of domain and range are the same, 
such a homomorphism is necessarily an isomorphism, i.e., 
we can find suitable bases in each cell and for each of 
the matrix algebras B q such that T\j(A) = UAU* for a 
unitary operator 



U: (g)C d 

xEU 



iC«(«) 



(16) 



Such a unitary operator by itself does not fix a QCA, 
because if we only take B q as an abstract matrix alge- 
bra, we still need to specify how r 1 _ q B q is contained in 
A(0 + 1) or, equivalently, to specify the isomorphism 
01 ® q T i-q®q with *4(D + !)■ This will be affected by 
another unitary operator 



V 



I c n(q) 



(17) 



ceD+i 



Any pair of unitaries (£7, V) according to Eqs. I|16I17|) 
specifies a transformation T on local algebras, which sat- 
isfies all requirements for a cellular automaton, except 
translation invariance: T only commutes with even trans- 
lations. The scheme in one and two lattice dimensions is 
visualized in Figures 17191 

□□□□□□□□□□□□ 

□□□□□□pnpngn 

'~v~ J v ~v~ J y ~v~^ ^~v^ "T" 
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FIG. 8: Generalized Margolus Scheme in s = 2 dimensions. 
The four subalgebras B q are symbolized by the shapes in Fia.Wl 
The tesselation is due to M.C. Escher (1956). 




FIG. 9: The four shapes representing the algebras B q in Fia.UH 
The arrows indicate the appropriate quadrant vectors q. 



FIG. 7: Generalized Margolus Scheme in s = 1 dimension. 
The algebras B+i and B-\ are symbolized by the different size 
cells in the intermediate step. 

If we want T, as constructed from unitaries U and V, 
to be a proper cellular automaton with full translation 
invariance, there will be additional conditions on these 
unitaries. Unfortunately, these conditions arc not easily 
written down and solved in the general case. We also 
note that U and V are not uniquely determined by the 
automaton: we have the freedom to choose a basis in 
every C™^. Changing this basis amounts to a cellwise 
rotation included in [7, which is immediately undone by 
the V-step. The key feature of automata in a partitioned 
scheme is structural reversibility, as described in the fol- 
lowing Lemma. 

Lemma 5 Let T be a homomorphism constructed from 
unitaries (U, V) in the generalized Margolus scheme. 
Then T is invertible, and T^ 1 is also a generalized Mar- 
golus automaton. Moreover, ifT commutes with all (not 
just even) translations, then both T and T -1 are nearest 
neighbor QCAs. 



Proof. The Margolus unitaries defining T -1 are (V* , U*). 
To check the localization properties, it is helpful not to 
think of these transformations in the apparently time- 
asymmetric scheme of Fig.|HI but to take B q as an algebra 
localized in the intersection of the cubes □ and (□ + q), 
i.e., as localized at the cell (1 + q)/2. Note that T _1 is a 
two-sided inverse, and hence uniquely determined by T. 

Therefore, if T commutes with all translations, so does 
T _1 . It remains to check that if T is a Margolus au- 
tomaton commuting with translations, it is actually a 
nearest neighbor automaton, i.e., a QCA with neighbor- 
hood scheme Af from 112(1. Since, for every x G □, we 
have £ (□ — x), we have 

T(Ao) C T(A(n - x)) C (g) A{D + q - x) . (18) 

<?eQ 

Since this is valid for all x £ □, we have T(Aq) C A(Af), 
with 

N= p| [j{U + q-x)=M (19) 
ie □ qeQ 
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□ 



B. Main Theorem 

Theorem 6 Let T be the global transition homomor- 
phism of a nearest neighbor quantum cellular automaton 
on the lattice If with single-cell algebra Aq — Aid- Then 
T can be represented in the generalized Margolus parti- 
tioning scheme, i.e., T restricts to an isomorphism 

T:A(n)^(g)B q , (20) 

where for each quadrant vector q G Q, the subalgebra 
B q C A(C1 + q) is a full matrix algebra, B q = A4 n ( q ). 
These algebras and the matrix dimensions n(q), which 
satisfy Eq. jl4j ), are uniquely determined by T . 

Combining this with Lemma we get 

Corollary 7 The inverse of a nearest neighbor QCA ex- 
ists, and is a nearest neighbor QCA. 

Note that this result is in stark contrast to the classical 
situation. In the classical case the inverse of an injective 
CA is a CA, i.e., locally invertible, but it is a highly non- 
trivial matter to determine the neighborhood scheme of 
the inverse, which can be much larger than the neigh- 
borhood of the automaton itself. This is not a contra- 
diction with the observation that every classical CA can 
be quantized (see Section fill B() : In order to construct a 
QCA from a classical CA, we needed the neighborhood 
size of the inverse. 

The proof of the Theorem will be given in subsec- 
tion llV Dl The key idea is to construct B q C A(0 + q) ex- 
plicitly from the inclusion T(A(D)) C (g> 9 A(D+q). This 
construction, which will also be useful independently, will 
be described in the next section. 



C. Support Algebras of Local Rules 

By definition, the transition rule To maps one cell alge- 
bra into a tensor product of neighboring ones. Therefore 
we need to investigate just how one subalgebra can sit 
inside a tensor product of others. 

Consider a subalgebra A C B\ ® B 2 of tensor product. 
For the moment let us forget about the multiplication 
laws, and just consider these as vector spaces, with the 
tensor product known from the (multi-) linear algebra of 
finite dimensional vector spaces. Then we can expand 
each a E A into a sum a = ^1 ® . But we might 

(i) 

get by just using a small subset of operators by . The 
smallest subspace of B\ sufficient for these expansions will 
be called the support of A on the first factor, and will be 
denoted by s(A,B\). For a more formal definition note 
that each a € A can be expanded uniquely in the form 



a = £^ a v ® e M> wnere i e ti} is fixed a basis of B 2 . Then 
s(A,B\) is the linear span of all a M in this expansion, 
and is clearly independent of the basis {e M } chosen for 

B 2 . Then it is clear that b^' £ s(A, Bi) indeed suffice to 
expand every a £ A, i.e., 

Acs(A,Bi)®s(A,B 2 ) cBi<g>B 2 . (21) 

The analogous relation for A contained in a tensor prod- 
uct of more factors is seen in the same way. 

Now we remember the algebraic structure: s(A, Bi) 
is a linear subspace of a C*-algebra, so we can define 
the support algebra S(A, Bi) of A C Bi on one factor 
Bi as the subalgebra of Bi generated by the elements of 
s(A, Bi) [24[. Then we also have 

A C S(A,Bx) ®S(A6 2 ) c B x ®B 2 . (22) 

Note that since any observable algebra A is closed under 
adjoints, so is S(A,Bi). The crucial fact we need about 
such inclusions is the following: 

Lemma 8 Let Ai C B\ ® Bi and A 2 C Bi ® B% be 
subalgebras such that A\ ® I3 and li C§) Ai commute in 
Bi®B 2 ®B^. Then S(Ai,B 2 ) and S(A 2 ,B 2 ) commute 
in B 2 . 



Proof. Pick bases {e M } C £>i and {e' v } C B 2 , and let 
a G A\ and a' £ A 2 . Then we may expand uniquely: 
a = e^^ti^ and a' = ^}2 V a' v ®e' v . Then by assumption 

= [a ® 1 3 , li ® a ']=Y^ e» ® a' v ] ® e l ■ 

fJ,U 

Now since the elements e^®e' v are a basis of B\®Bz, this 
expansion is unique, so we must have [a^a'J = for all 
/i, v. Clearly, this property also transfers to the algebras 
generated by the and a' v , i.e., to the support algebras 
noted in the Lemma. □ 



D. Proof of the Main Theorem 

We apply the construction of support algebras to the 
inclusion 

T{A(n))c®A(D + q) (23) 
<? 

and define 

B q = 8(T{A(p)),A(p + q)). (24) 

As a finite dimensional C*-algebra, each B q is isomor- 
phic to B q = © A1 A4 n ( g M ) (See Proposition H*T|l . Now 
T(A(0)) is homomorphically embedded into B q , with 
T(l) = 1. Hence by Proposition^|we know that for any 
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choice of summands fi q we must have that d 2 , the ma- 
trix dimension of A(\3), divides Y\ q n{q, n q ). This gives 
a lower bound on the block sizes. 

In order to get an upper bound, consider the support 
algebras contained in some shifted cube, such as □ + 1. 
These are 

s(r(A(D + l-q)) , A(D + 1)) 

= T 1 _,s(r(^(n)),.4(n + ? )) 

= r!_,(B,) (25) 

Since the A(0 + 1 — q) commute, so do their images 
under T and, by Lemma|Sl so do the algebras r 1 ^ q (B q ) C 
A(0 + 1). However, we do not know a priori that these 
algebras arc contained in „4(D + 1) like a tensor product: 



When 



G B q is a central projection onto one of the 



q ^q^ q 



matrix blocks of B q , it is clear that the product \[ 
is a central element of the algebra generated by the B q , 
but it might be zero. On the other hand, there must be 
some choice of blocks fi q , for which this is non-zero, and 
for this combination Y\ q M n ( q ,^ q )i which is isomorphic 
to A4 n with n — 71(5, , is a direct summand of 
Y\ q £> g C A(D + 1). Hence we get the inequality 



d 2 " > Y[n(q,fi q ) 



(26) 



On the other hand, by the first step, the left hand 
side divides the right hand side of this inequality, so the 
we must have equality. This also implies that only one 
summand can be present in B q , so we get B q = M n ( q ) 
with n(q) = n(q,fi q ). That T is an isomorphism from 
A(0) onto @ q B q follows by a direct dimension count, 
since a *-homomorphism between full matrix algebras of 
equal dimension can only be zero or an isomorphism. □ 



E. QCAs with abelian neighborhood 

In a sense, Theorem [5] gives a complete constructive 
procedure for QCAs in terms of the two unitary oper- 
ators U , V with a constraint. Unfortunately, however, 
it does not seem to be easy to give a general solution 
of the constraint equations expressing the translation 
invariance (rather than the invariance by even transla- 
tions). Therefore it is suggestive to repeat the analysis 
of support algebras also on the single cell level. Setting 



we have 



V x = S(T (A ),A X ) 



T {Ao) C (g) V x 



(27) 



(28) 



It is clear that since To(Aq) is non-abelian, at least one of 
the algebras T> x must also be non-abelian. The simplest 
case in this regard will be when all T> x are abelian and 



commute with each other, except one, say T>q, which then 
has to isomorphic to the full cell algebra Aq by Propll2l 
Since all T> x commute, we can jointly diagonalize them 
and this fixes a canonical basis for every cell. When 
denotes the basis vectors, we can write the local transi- 
tion rule as 



T (A) = Y,U( m )*AU( W ) 



lMx)(Mx| , (29) 



where the sum runs over all tuples fi^f of basis labels 
\x x for x £ Af , x ^ and, for each such tuple, U (///v) 
is a unitary operator. Thus Tq describes a conditional 
unitary operation on cell 0, where the conditioning is in 
some fixed "computational basis" . 

We now need to analyze the constraints on these uni- 
taries needed to make this homomorphism To a local 
transition rule. As a first step we look at the simplest 
case: 

Lemma 9 Let U^, V v € A4d be unitary operators (ji, v — 
1, . . . , d) such that, for all A, B £ Md, 



J2[u;au^\^)(^\, \ V ){v\®v:bv v 



= 



(30) 



Then there are unitary U, V £ Add such that, for all \i, 
U*U^ and V*Vfj, are diagonal. 

Proof. Let us abbreviate A^ = U^AU^, B u = V*BV Vl 
and take the matrix element of equation (|30|l in the prod- 
uct basis vectors {af3\ ■ ■ ■ \^5). This gives 

(a\A p \i) (/3|B 7 |J) = (a|A 5 | 7 ) {(3\B a \5) 

Now set B = V a \(3')(5'\V*, with [3' j= /3. Then B a = 
\[3'){8'\, and {(3\B a \5) — 0, and the right hand side van- 
ishes. Hence, for every A and every 5' 

(a\A p \j) (P\V;V a \f3')(8'\V:V^\6) = 0. 

Now the first factor can be made non-zero by an appro- 
priate choice of A, and the third factor can be made non- 
zero by choosing 5', because the unitary operator V*V^ 
cannot annihilate \6). It follows that ((3\V*V a \(3') = 
vanishes for all indices, or V*V a is diagonal for all a,j. 
Hence the Lemma follows with V = V\, and the state- 
ment for U follows by symmetry. □ 
Let us apply this Lemma to the one-dimensional near- 
est neighbor case. Then a local rule of the form (|29|l can 
be written as 



T (A) = ]T |/x)(Ml «8 U^AU^ ® \u)(u\ 



(31) 



The remaining commutation condition is, for arbitrary 
one-site observables A,B, 

= [T (A)®1,1®T (B)} 
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Hence we can apply the Lemma to the commutator sep- 
arately for every pair ji, u' , and we find that up to a 
common unitary all have to commute. Again, up 
to a cell-wise unitary rotation, we can choose the com- 
mon eigenbasis of the as the same basis in which the 
conditions are written, i.e., 



U, 



(33) 



with some phase function u depending on three neigh- 
boring basis labels. However, this phase function u is not 
arbitrary: unitaries of the form l|33l) in this way may still 
fail to satisfy (|3~2"j) . If we insert A = |a)(6|and.B = \a')(b'\ 
we get the functional equation 



u(fiba') u{Wv') _ u(fibb') uiab'v') 



(34) 



u(fiaa') u(ba'v') u(fiab') u(aa'v') 

Since we want to classify solutions up to a cell-wise ro- 
tation, we can take one of the unitaries U^ v to be the 
identity, say U\\ = 1, or u(lxl) = 1. Moreover, an over- 
all phase of is irrelevant, and we can choose this so 
u{ixlv) = 1. Then in l|34l) we take a = V = v 1 = 1, which 
gives 



u(fiba') — u(fibl)u(ba 1) 



(35) 



Thus u is already determined by the two-variable func- 
tion (a, V) *—> u(a, b, 1). It is easy to check that any choice 
of this function yields a solution of l|34|) via (|35|l . 
We can summarize the result as follows: 

Proposition 10 Let T be the local transition rule of a 
QCA such thatT>\ andT>-\ are both abelian. Then, with 
respect to a basis in which these algebras are diagonal, 
there is a phase gate on C d ® C^: 



U = ^2u(a,b)\ab)(ab\ 



(36) 



ab 



normalized such that u(l, b) = u(b, 1) = 1, and a one-site 
unitary V such that 

T (A) = X*(t ® V*AV ® 1)X with 
X = (l7®ia)(li®C0 . 

F. Unilateral Automata 

Another case of automata in one dimension, which 
can be characterized completely, are automata with 
neighborhood scheme J\f = {0,1} (or, symmetrically, 
M = { — 1,0}). The analysis is almost identical to 
that of Theorem El and gives full matrix algebras T>i — 
S(Tq(Ao), Ai) = Aim and such that hqUi = d. As in 
the case of the Theorem, the local rules combines an 
arbitrary unitary U : C d — > C™° ® C™ 1 with a unitary 
V : C" 1 ®C"° -» C d . We only mention these to state the 
following classification of the simplest case: 



G. Nearest neighbor qubit automata in one 
dimension 



Consider the right and left support algebras 
T>-i, T>q, T> + i as in (I27|l . These are subalgebras of the 
2 x 2-matrices, which leaves three possibilities: each of 
these algebras can either be trivial (£>.; = CI), an abelian 
two-state algebra (isomorphic to the diagonal matrices, 
or the full algebra 7^2- 

Suppose that at least one of the algebras 2?±i, say 
is trivial. Then we have a unilateral automaton. Since 
n^ni = 2 we must have either ng = 2, n\ = 1, a cell-wise 
unitary rotation or uq = 1, n\ — 2 a right shift, possibly 
combined with a cell-wise rotation. 

Suppose that none of the algebras T>±\ is trivial. Then 
because X>_i commutes with T>+\ neither algebra can be 
the full matrix algebra, since that would force the other 
to be trivial. It follows that both are abelian, and com- 
mute. Hence after a basis change (by another cell-wise 
rotation) we can take 2?_i = T> + i as the algebra of di- 
agonal 2 x 2-matrices. This brings us into the situation 
of Section IIV El and we find a phase rotation. Choos- 
ing the normalization in Proposition 1101 we have only 
one free parameter left, i.e., we have an automaton built 
from commuting unitaries U x (cf. Section [illfl . which are 
phase gates 



U x 
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The classification is hence complete. It is represented 
in FigUni 
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FIG. 10: All nearest neighbor qubit automata arise by combin- 
ing cell-wise unitary rotations (C) with right shifts (R), left 
shifts (L), or phase gates (P). 

It is interesting to compare this with the classical case, 
which can be stated very similarly: then the only cell- 
wise operations are identity and global flip. Of course, 
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the possibility of phase gates does not arise, leaving 6 
classical possibilities. 

Beyond qubits, a good classification exists for Clifford 
automata in prime dimension p^ . However, for gen- 
eral three-level systems the classification is likely to be 
complicated, since there is already a host of reversible 
classical nearest neighbor automata. 

V. OTHER APPROACHES 

In this section we take a look at the various propos- 
als for defining QCAs, and show how they relate to the 
approach taken in this paper. 

A. Feynman, and Transition Quasi-Probabilities 

The idea of a quantum cellular automaton is clearly 
present in Feynmans's famous 1981 lecture 0. Although 
he suggests not only that a theory might and should be 
developed, and that it might even be taken seriously as 
fundamental physical theory, he does not actually de- 
velop a notion of QCAs in this article. The context in 
which he does write down a transition rule for a QCA, 
is where he discusses the possibility of simulating a QCA 
with a probabilistic cellular automaton. He emphasizes 
that this would involve something like negative transition 
probabilities and closes this section saying "... I wanted 
to explain that if I try my best to make the equations 
look as near as possible to what would be imitable by a 
classical probabilistic computer, I get into trouble" . 

Feynman's idea for making a QCA look as classical as 
possible is to replace "transition probabilities" by Wigner 
function-like "transition quasi-probabilities" . This ap- 
proach is an interesting contribution to the definition of 
irreversible QCAs, which is still plagued with problems. 
However, as Feynman is clearly aware, it fails, and it is 
instructive to analyze this failure in the case of reversible 
automata. 

The transition function of a classical probabilistic cel- 
lular automaton (for a set S of single-cell states) is a 
set of probabilities M(s'\sm) specifying the probability 
of finding a state s' € S, when the configuration of the 
neighbors at the previous time step is sjj. Using the 
"simultaneous and independent update rule" the prob- 
ability for finding a configuration s\ in a finite region 
A C Z s becomes 

M(s' A \s) = l[M(s' x \s A r +x ). (38) 

xeA 

Note that this is readily read as a statement in the 
Heisenberg picture: the probability for s' A is expressed 
as the expectation of a random variable in the previous 
time step, namely the right hand side of ijSSJ, consid- 
ered as a function of the variables s x . In the quantum 
case, M(s'\sj\f) would then become an observable in A.j\f, 
depending in a linear (and completely positive) way on 



an observable s' £ Aq. This is precisely a description 
of the local transition rule To. Note that the commuta- 
tion condition for local rules (Lemma[2l implies that the 
product is well-defined, independently of the ordering of 
the factors. In a more general quantum context, such 
as irreversible QCA evolutions for which the global rule 
will not respect the product, we cannot be sure of this 
property: additional information about the ordering of 
factors in l|38(l would have to be supplied, but even an 
ordering fixed by some convention would not prevent the 
evolution from sometimes taking hermitian elements to 
non-hermitian elements. 

This latter problem: the ordering of factors and the 
hermiticity is neatly solved by Feynman's approach of 
quasi-probabilities. He expands all qubit operators in 
a special basis of four hermitian operators F(£) <E M2, 
{++,+--+,— }: 

F -( (1 + U)/2 "p-to)/*] (uv-±l) (39) 
* uv -\v\l + iu)/A (l-u)/2 ) ,lM, "" ±1J ' 

The expectations / p (£) = tr(pF(£)) of thes e op erators 
are the analogs of the Wigner function (see [25|, |2g for 
later elaborations on Wigner functions in finite dimen- 
sion). Of course, by taking tensor products of these op- 
erators we get Wigner functions for multi-qubit systems. 
Now we can expand the transition rule for a single cell 
in Wigner operators: 

To(F(vo)) = ]T Mfa,|&v) (g) F(£ x ) , (40) 
Sm xeAT 

with the transition quasi-probabilities M(?7o|Cv)- These 
are usually not positive, but this is a minor inconvenience 
as long as the physical transition operator T is positive. 
Equation (|40l) is just an expansion of the local rule in a 
basis of Wigner operators, so a local rule is completely 
equivalent to a set of transition quasi-probabilities. The 
difference between the quasi-probability approach and 
ours is how the global rule is constructed. Let us con- 
sider, for simplicity, the image of a two-site observable 
F(j]i) (8)^(172) under a one-dimensional nearest neighbor 
automaton. According to i|38|l . transition probabilities 
must be combined as 

Tquasi 

(F( m ) ® F( V2 )) 
= Yl l£o,£i,6)M(77 2 (41) 

^F{^)®F(^)®F(^)®F{^) . 

We can extend this to arbitrarily large configurations to 
get the time evolution of an automaton in the Heisen- 
berg picture. Since the transition quasi-probabilities are 
real, this evolution will automatically preserve hermitic- 
ity, and since F(£) = 1 all normalization and local- 
ity properties will automatically come out correctly. Of 
course, there will be no operator ordering problem, since 
we multiply at the level of functions, and all this will 
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work for probabilistic irreversible and reversible transi- 
tion rules alike. 

On the other hand, under the homomorphism T the 
tensor product F(r?i) <g> F(?7 2 ) is evolved to 

T(F(rn) 8> Ffa)) 
= Yl M (Vi\^,^,^)M(r, 2 \^,^,Q (42) 

€o>£i.£2>€3 

x F(£o) ® F(fc)F(£) ® F(&)F(&) ® F(^) 

Then gTJ=g21>, iff F(£)F(r)) = 6 0l F(£), which is ob- 
viously true for the minimal projections in a classical 
observable algebra, but obviously false for the Wigner 
operators. So the automata studied in this paper are not 
(or not in general) representable as quasi-probabilistic 
CAs. 

Feynman uses the approach based on ((41(1 only as a 
caricature of a quantum process. He points out that neg- 
ative transition quasi-probabilities exclude a simulation 
of the global time step 1(38(1 by successive independent 
random trials. For him this indicates the increased com- 
plexity of quantum computation. 

Whether or not this approach works as a definition 
of (possibly also irreversible) QCAs hinges on the ques- 
tion of positivity. Non-positive transition probabilities 
would not be serious if they lead to the construction of 
a legitimate (i.e., completely positive) global transition 
rule. Unfortunately, however, there is no indication that 
quantum positivity improves in the passage from local 
to global rule. This can be seen already for a simple 
phase gate array, i.e., (P) in Fig. ^] For generic phase 
angle (p neither the local transition nor the global tran- 
sition e.g., the expression 1)42(1 . have positive transition 
quasi-probabilities, although, of course, they correspond 
to completely positive operations. On the other hand, 
the two-site transition rule (f 4 If) takes some positive op- 
erators into non-positive ones. 

An interesting special case occurs for phase angle 
(j) = 7T. In that case, the local rule does give rise 
to positive, and even deterministic transition quasi- 
probabilities. They belong to a classical deterministic 
CA, which like the phase gate QCA is its own inverse. 
But if it is applied as in 1)41(1 , it violates positivity. 

B. Watrous et al. 

One of the first serious attempts at the definition of 
QCAs was by Watrous It is based on another "quan- 
tization" of the transition probability formula 1(38(1 . in- 
spired by Feynman's notion that in quantum theory one 
must replace probabilities by amplitudes. Thus one tries 
a product formula like ((38(1 for the transition amplitudes, 
presumably defining in this way the unitary transition 
operator for the whole process: 

U(a\b)=H u(a x \bj>s +x ) , (43) 

X 



where a, b are classical configurations, labelling the basis 
states of the QCA, and b^ +x is the configuration b re- 
stricted to the neighborhood of x. The function u plays 
the role of the local transition rule. Basically the same 
definition is also used in van Dam |2j , where it is phrased 
as an assignment of a a product vector to every basis 
state in the computational basis. Further work in this 
approach is to be found in |9(, and in the textbook |8|], 
or a recent introduction |27f . For the sake of discussion 
let us call an automaton defined by ((43(1 a WQCA. There 
are several problems with this formula: 

1. The infinite product may not be defined. This 
may be resolved by either introducing a "quies- 
cent state" which is invariant under the evolution, 
and considering only superpositions of configura- 
tions which are quiescent outside a finite region 
0, EljOr to look at periodic boundary conditions 
only |3 • Either approach works, but none of these 
workarounds is necessary in our definition. 

2. U from ((43(1 has no reason to be unitary, and most 
of the time it isn't. In other words, there is no 
straightforward way of characterizing those local 
transition amplitudes u for which the formula does 
indeed define an isometric (or stronger: a unitary) 
operator, in which case the rule is called "well- 
formed" (or unitary). Whereas positivity and nor- 
malization of the local rule M on the right hand 
side of l(38|) guarantee the corresponding properties 
for the global evolution, no equally simple criterion 
exists for well-formedness (but see [2j for algorithms 
in the one-dimensional case). 

3. Many unitary operators are not of the form 143(1 . 
To begin with, the definition depends on the choice 
of a preferred basis. In general, the product of two 
WQCA unitaries need not be a WQCA (with larger 
neighborhood), and the inverse of a WQCA may 
fail to be a WQCA. This makes it hard to build a 
general theory on this definition. 

To get an example of these phenomena note that a 
necessary condition for a unitary of Watrous form 
is that the computational basis states are mapped 
to product states. Moreover, if we follow the uni- 
tary operator by a site-wise unitary rotation, or 
precede it by a product of phase gates (i.e., intro- 
ducing additional phase factors independent of the 
output labels a x ), we stay in this class. But now 
consider a sitewise Hadamard rotation, followed by 
a product of phase gates. It is easy to verify that 
such a map does not take the computational basis 
states to product states (or Briegel's one-way com- 
puters would not work.) So this is not a WQCA, 
although it is the product of two WQCAs, and its 
inverse is also a WQCA . 

4. Even in the cases where the formula does work, 
such as, e.g., for the multiplication by local phases, 
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the size of the neighborhood in (|43|) is not the neigh- 
borhood scheme describing the propagation of ob- 
servable effects (the TV of our definition) , but rather 
the Af from Section IIII Al In fact, it is not even 
clear whether a WQCA is necessarily local in our 
sense: We did not manage to exclude the possibility 
that a local measurement after one time step might 
allow inferences about arbitrarily distant modifica- 
tions of the state. 

5. in the original paper |||, Watrous also looks at par- 
titioned automata. Since we prove that all QCAs 
can be obtained by a partitioning scheme, it might 
seem that WQCAdQCA, in seeming contradiction 
to the Example in 3 above. However, Watrous uses 
only a special case of partitioning, namely a uni- 
tary followed by a permutation of subcells. Only 
if we extend the class of WQCAs and include all 
their products, we get all QCAs. 

To summarize: Judged by the criteria in the introduc- 
tion, the notion of Watrous is not a satisfactory formal- 
ization of the idea of quantum cellular automata. What is 
lacking is the easy passage from local to global rules, and 
from global axiomatic to local constructive description. 
For all problems involving the iteration of automata, the 
fact that the class is not closed under composition and 
inverse, puts a premature end to studies based on WQ- 
CAs. 



C. Richter and Werner 

As mentioned earlier, the observable-based approach 
underlying our definition was first used in by one 
of us, but with a focus on the irreversible case. In that 
paper partitioning (dissipative cell evolution combined 
with permutation of subcells) was used to allow a free 
construction satisfying a global locality condition. 

It is not so clear what should replace the axiomatic 
definition in the irreversible case. A possible condition 
is to just replace the local rule To by a completely pos- 
itive map, and insist on commutativity as before. This 
does define a global evolution step 0, |28| . In view of 
the reversible case this would seem like a good definition. 
However, it is again not clear whether the composition of 
two such transformations will again be of the same kind. 
Nor is it clear how to describe the class of transformations 
obtained by several such "simultaneous independent up- 
date" steps. 

It turns out that the analysis of Theorem [5| can partly 
be repeated in the irreversible case, thereby reducing the 
possibilities somewhat. This line will be pursued else- 
where. 



D. Wolfram 

In a recent thick book [2^| S. Wolfram has argued 
that the that the universe might be a big cellular au- 
tomaton following simple rules (see also 30]). Wolfram 
uses only classical structures, expressing the belief, how- 
ever, that quantum structures might emerge from classi- 
cal rules generating sufficient complexity. Since Einstein 
failed with a program like that (he thought of overde- 
termined non-linear field equations, rather than CAs) it 
would be nice to see the details worked out. 



E. Quantum random walks 

The term "Quantum cellular automaton" has some- 
times been used |3lL l32t l33t for a unitary evolution 
of a particle on a discretized space. The total Hilbcrt 
space of such a system is £ 2 (Z S ), the space of square 
summable complex functions on the lattice, i.e., the di- 
rect sum rather than the tensor product of the one-cell 
spaces. The classical analogue of such a system is a single 
classical particle moving on a lattice, e.g., in a random 
walk. Therefore much better terminology to call these 
quantum systems quantum random walks, rather than 
cellular automata. Such systems have been proposed 
for purposes of quantum computation, in particular for 
search problems on graphs (see |35j for a review). 

Quantum random walks require localization properties 
not unlike those of QCAs. To see the connection it is 
interesting to consider the connection between classical 
random walks and CAs: A random walk can be seen as 
a special CA, with each cell either empty or occupied, 
started in a configuration with exactly one occupied cell. 
Of course, the overall CA dynamics should respect this 
constraint. Because the CA rule is local we can then 
also put arbitrarily many particles on the lattice, and the 
dynamics is well-defined by the random walk, as long as 
the particles do not collide. What happens on collision is 
a piece of information which must be supplied in order to 
make a random walk into a CA, i.e., in order to pass from 
a random walk to a (possibly "interacting" ) diffusion. 

Consider a QCA with a special "empty" state speci- 
fied for the single cell algebra. This means that we can 
define a global quantity particle number, which ought to 
be conserved by the QCA. Technically this is the infinite 
sum of O's and l's, and not a well defined observable. 
However, we can consider this formal sum as a lattice in- 
teraction generating the time evolution of "gauge trans- 
formations" , acting as an infinite product of unitaries, 
as in Section IIII Al For a QCA commuting with such 
transformations, the "one-particle" Hilbert space can be 
defined, and the QCA dynamics restricted to this sub- 
space is a unitary evolution of the random walk type. 
Note that we have not assumed that the dimension d of 
the one-site algebra is 2, i.e., occupied cells ("particles") 
may have an internal structure. This turns out to be 
necessary: a standard classical random walk, with only 
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the empty /occupied distinction and nearest neighbor in- 
teraction cannot be reversible. Similarly, a unitary on 
£ 2 (Z S ) cannot be strictly local [jy]], and we do need in- 
ternal states [3^. If there is only one particle, this is the 
same as saying that we have only the empty/occupied 
distinction for the cells, but we have also a "quantum 
coin" which helps determining the steps. 

The standard model of a quantum random walk |3fij 
uses a qubit coin, i.e., three states for the QCA. All ran- 
dom walks with such a coin are parameterized by a sin- 
gle unitary 2 x 2-matrix U: we can describe the internal 
states ("chirality" 0) as "go right" and "go left". A 
trivial, but globally well-defined evolution step is defined 
by following these instructions. The general case arises 
by following this with a sitewise unitary rotation by U, 
representing the quantum coin flip. 

Now the question arises: can we consider such random 
walks as the one-particle component of a QCA allowing 
arbitrarily many particles? Indeed this is possible, even 
in many ways, and it would be very interesting to classify 
all possibilities, hence all "interactions" . One possibility 
is to second quantize the random walk, which leads to a 
Boson system allowing, at each site, an arbitrary number 
of particles. It is clear from the formalism of second quan- 
tization that all the locality properties required of a QCA 
are then satisfied. Clearly, this is the non-interacting op- 
tion. If we insist on finitely many states per cell we in- 
troduce some kind of hard core interaction. A trivial way 
of doing it with four states is this: Consider two infinite 
qubit chains, called the left moving and the right moving 
chain. The "free" time evolution S is indeed shifting the 
two chains separately according to this description. Now 
at each lattice site we take the right moving and the left 
moving one-site algebras together to define a single cell of 
the QCA, with the four states labelled 'empty', 'R', 'L', 
and 'RL'. We then select a unitary U leaving the empty 
and the doubly occupied state invariant, but shuffling 
the 'R' and 'L' states as before. Clearly, the combination 
of the free evolution S with the sitewise application of 
U defines a QCA whose one-particle sector is the given 
random walk. 



F. Spacetime localized algebras 

A quantum field theory can be considered from two 
equivalent points of view. On the one hand one can con- 
sider the fields at some time fixed time as the basic dy- 
namical variables, e.g., the Cauchy data at time zero. On 
the other hand, for the discussion of relativistic causal- 
ity, it is convenient to consider as fundamental the fam- 
ily of algebras 21(0) associated with the measurements 
in some space-time region 13| . Similarly, one can look 
at a QCA from these two points of view, the "Cauchy 
data" point of view being what we described so far. For 
the space-time view consider the extended lattice Z s+1 
of space-time points (x,t), with x € Z s , the same spatial 
lattice as before, and t £ Z a discrete time point. 



The global C*-algebra will be the same as before, but 
we introduce additional subalgebras, namely 

A x , t = T-\A x ) . (44) 

Since T commutes with lattice translations we can de- 
fine a joint set of space-time translations T Xjt , combin- 
ing a lattice translation by x with T~ l . For any subset 

c Z s+1 we define 21(0) as the C*-algebra generated 
by all the A x ,t with (x, t) £ 0. In order to study the lo- 
calization properties of of these algebras, let us say that 
a sequence of lattice points 

(x ,to), (xi,t + l), (x 2 ,t + 2),...,(x N ,t N ) (45) 

is (forward) timelike, if Xk+i — Xfc € M for all k. Then we 
define two regions 0i,02 C Z s+1 to be spacelike sep- 
arated (notation: 0i)(02), if no timelike curve passes 
through some point of 0\ and also through some point 
of 02. Moreover, we say that 0\ is causally dependent 
on 02 (notation: 0\ < 02), if every timelike curve which 
passes through (any point of) 0\ also passes through 02 . 

Then it is clear that 0iX02 implies that 2l(0i) and 
21(02) commute elementwise. Moreover, if 02 is a finite 
set of lattice points (x,t), all with the same t, and if 

01 < 2 , then 2l(0i) C2t(0 2 ). 

Conversely, if we have a net of local algebras 21(0), 
defined for arbitrary finite subsets of the spacetime- 
lattice, if the spacetime translations act by automor- 
phism of the total algebra such that T( x jt )(2l(0)) = 
21(0+ (x, t)), and that the above locality properties hold, 
then the time zero algebras form a QCA in the sense of 
Section |nl 

This way of viewing QCAs may indeed be useful for 
making the connection to relativistic field theories, or for 
a detailed study of the growth of localization regions. 

G. Non-commuting cells 

As is standard in quantum theory we described the 
cells of the automaton as subsystems in the usual ten- 
sor product sense. In a slightly more axiomatic style 
we could have postulated, that arbitrary measurements 
in separate cells can be carried out jointly, implying the 
commutation between the different cells. While this is 
certainly very natural, more general commutation rules 
between cells may be considered. An obvious example are 
anti-commutation rules, i.e., we could think of a Fermi 
gas on a lattice. This would, of course, change the re- 
quirements on local rules, but it is quite clear how to 
adapt our definition to that case. 

Another very interes ting deviation from commut- 
ing cells is studied in |37|. in connection with non- 
commutative analogs of 2-surfaces with constant negative 
curvature, the sine-Gordon equation, and 2-dimcnsional 
lattice systems in magnetic field (Hofstadter butterfly). 
Roughly speaking the structure investigated has has a 
single variable, say a unitary Q Xl at each cell Non- 
commutativity comes in, because any pair of neighboring 
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unitaries forms a discrete Weyl system. At distances > 2 
the variables commute. It is clear that (reversible) dy- 
namical rules can be described exactly along the lines 
of our definition, as automorphisms respecting this alge- 
braic structure, and also the localization (up to a finite 
enlargement of localization regions). Of course, this is 
not the place to explore such structures, and we point 
again to the book |37j for more material. 



APPENDIX A: FINITE DIMENSIONAL 
C*- ALGEBRAS 

Almost all the hard work in any textbook on C*- 
algebras goes into aspects of the theory, which become 
entirely trivial for finite dimensional C*-algebras. Of 
course, some algebras in this paper are infinite dimen- 
sional, notably the quasi-local C*-algebra describing the 
infinite system. But the key arguments use only the fi- 
nite dimensional structure. Therefore we will give in this 
Appendix a quick summary of C*-algebra theory as it 
applies to the finite dimensional case. 

In order to make the paper more accessible to commu- 
nities in which algebraic terminology is less current (e.g., 
most theoretical physicists, and classical computer scien- 
tists) we start out with an extended glossary, in which 
the basic notions are defined and some basic facts are 
noted. This will be followed by some structure theorems 
which we need in the body of the paper. Of course, all 
this cannot replace a serious textbook. We recommend 

Emm. 



1. C*-Glossary 

The operations making up the abstract structure of 
C*-algebras are inspired by those known from algebras 
of operators on a Hilbert space. In fact, every algebra of 
operators on a finite dimensional Hilbert space, which is 
also closed under taking adjoints, satisfies the definition, 
and every abstract C*-algebra is isomorphic to an op- 
erator algebra. The following list of relevant operations 
and concepts may serve as a glossary. The algebra under 
consideration is usually denoted by A. 

• Addition and multiplication by complex scalars. 
This makes A a vector space over C, which we as- 
sume to be finite dimensional from now on. 

• Multiplication. We denote the product by AB G A, 
when A, B G A. The product is distributive and 
associative (but not necessarily commutative). If 
the algebra is commutative or "abelian" the system 
under consideration is classical. 

• The adjoint or "star operation" denoted by A* G A, 
when A £ A. This is conjugate linear (or "an- 
tilinear"), which means that (XA)* = XA* , and 



{A + B)* = A* + B* . The adjoint satisfies (AB)* = 
B* A* . Physicists often write A^ for the adjoint. 

• The norm, which is a positive number \\A\\ asso- 
ciated with each A € A. With respect to the al- 
gebraic structures, the norm satisfies + < 
\\A\\ + \\Bl \\XA\\ = |A| ||A||, \\AB\\ < \\A\\\\B\\ 
and ||AM|| = ||A|| 2 . ||A|| = implies A = 0. In 
contrast to the general case, the norm is uniquely 
determined by the algebraic structures. 

• An identity. In contrast to the general case, in finite 
dimensional C*-algebras there is always a unique 
element II satisfying 1A = At — A. 

• A homomorphism between C*-algebras is a map 
$ : A — > B between C*-algebras, preserving the al- 
gebraic structures. That is, <E> is linear, Q(AB) = 
$(A)$(B), and $(A*) = $(A)*. The latter prop- 
erty is sometimes emphasized by speaking of *- 
homomorphims. Homomorphisms $ : A — > A are 
called endomorphisms of A, and if an inverse ho- 
momorphism exists, $ is called an isomorphism or 
an automorphism (if A = B). For general homo- 
morphisms, $(1) is a projection in B. If $(1) = 1, 
we call it unital. 

• An ordering. We write A > 0, if there is a B 6 A 
such that A — B*B. The set of positive elements 
is a convex cone in the set of hermitian (A* = A) 
elements. In an operator algebra, the positive el- 
ements are precisely those hermitian ones with all 
eigenvalues non-negative. 

• The center of A is the subalgebra Z(A) C A of 
elements Z such that ZA = AZ for all A G A. 

• A state on A is a linear functional u> : A — > C 
which is positive (i.e., A > implies uj(A) > 0) 
and normalized (i.e., w(l) = 1). 

• A trace on A is a positive linear functional r such 
that t(AB) = t(BA). 

• a positive linear functional u> is called faithful if 
A > and w(A) = imply A = 0. Every finite 
dimensional C*-algebra has a faithful trace. On 
an operator algebra, the usual trace of operators 
(which we will denote by tr) is an example. 

• Given a faithful trace r, every positive linear func- 
tional can be written as u>(A) = r(pA), for a unique 
p G A, p > 0, which is called the density operator of 
u> with respect to r. oj is also a trace iff p G Z(A). 

• An element P G A is a called a projection if P* = 
P = P 2 . It is called a minimal projection if for any 
projection Q, Q < P implies Q = or Q = P. 

• The direct sum A = ® A^ of a finite collection 
of finite dimensional C*-algebras A^ is the vector 
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space direct sum, i.e., elements are tuples of com- 
ponents G Afj,, with componentwise algebraic 
operations. In operator algebras each term in the 
sum corresponds to a diagonal block in a block ma- 
trix decomposition. 

• The tensor product A = (^) A^ is the vector 
space tensor product, with product and adjoint de- 
fined as the unique linear (resp. conjugate linear) 
extensions of (^^((g^-B^) = ® M (A^S^) and 
= ^ n °P era t OT algebras one forms 

this product by taking first the tensor products of 
the underlying Hilbert spaces, and taking the alge- 
bra generated by all tensor product operators. 

The first four items on this list are the definition of C*- 
algebras. Note that order and unit are explicitly defined 
in terms of the algebraic structure (linear operations, 
multiplication and adjoint), and the norm is also defined 
explicitly as 

||;4|| - mf{A > 0\3 B A* A + B*B = A 2 1} (Al) 

It might thus seem superfluous to list the norm among the 
defining elements. In the infinite dimensional case it is 
needed, of course, to formulate the topological complete- 
ness requirement (and completeness is needed in turn to 
construct B in <|Aljl L However, even in the finite dimen- 
sional case the implication (\\A\\ = 0) {A = 0) car- 
ries non-trivial information by excluding the existence of 
nonzero elements Ai such that A$ Ai = 0. 

2. C*-structure 

Consider a single Hermitian element A G A. Then 
since A is finite dimensional, the powers A n must be lin- 
early dependent, i.e., there is a characteristic polynomial 
P(A) — Ylk c kA k — 0. From this one readily constructs 
polynomials pi such that pe(A) is a projection, and 

A = Y,aiPi(A) (A2) 

l 

where ai are the distinct roots of p(a) = 0. This is called 
the spectral theorem. It implies, in particular, that any fi- 
nite dimensional C*-algebra has many projections (which 
may fail in infinite dimension). 

This fact will be used in the following fundamental 
structure theorem. Recall that by Ai n we denote the 
algebra of complex n x n matrices. 

Proposition 11 Every finite dimensional C*-algebra A 
is characterized uniquely up to isomorphism by a finite 
sequence tlx > n% > • • • > 1 of numbers such that 

A = ($M n ^. (A3) 



The basic idea of the proof is to consider the center 
of A, which is a finite dimensional abelian C*-algebra. 
The minimal projections z^ of the center decompose the 
algebra into a direct sum A = z^A, in which each 
of the summands has trivial center. The building blocks 
z^A are then seen to be isomorphic to full matrix alge- 
bras .M njj , where is the maximal number of mutually 
orthogonal projections in z^A. 

Proposition 12 If $ : Aid — > M M n (n) is a *- 
homomorphism such that <&(!) = 1, each n(fi) has to 
be divisible by d. 

By considering the composition of the given homomor- 
phism with the projection onto one summand, which is 
also a homomorphism, we can consider the case of a sin- 
gle summand. Thus $ becomes a representation of Aid 
on a Hilbert space of dimension n(/x), which can be de- 
composed into irreducible representations. It is a basic 
property of Aid, however, that all its irreducible repre- 
sentations are unitarily equivalent to the defining repre- 
sentation on C d , so n{n) must be d times the multiplic- 
ity (number of isomorphic irreducible summands) of this 
representation. 

Proposition 13 If A S M % and B = ©„ M m{v) 
are commuting subalgebras of B(7i), the algebra AB is 
also decomposed into direct summands, each of which 
arises by multiplying a summand form each of the al- 
gebras. These occur with integer multiplicities r^„ > 
such that 

'^^ / r llu n(fi)m,(b') — dimH . (A4) 

Let A^ G M n ^) G A and B u G M m ( v ) G B be ele- 
ments from the respective blocks. Then A^B^ i— > A^B V 
is a representation of Ai n ^). m [ v ) = Ai n ^ (S> -M OT („) on 
Ti., which may however be zero since we cannot guarantee 
that it preserves the identity. The rest of the argument 
is as for the previous proposition. 
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